EPJ manuscript No. 

(will be inserted by the editor) 



OS 

o 
o 

(N 

O 
D 

Q 

(N 



CZ5 

O 

> 

q 
o 

>V 

^ : 

Oh. 



> 

(N 
(N 

(N 

On 
O 



Diffusion and Localization of Cold Atoms in 3D Optical Speckle 

A. Yedjour 1,2 and B.A. van Tiggelen 1 

1 Laboratoire de Physique et de Modelisation des Milieux Condenses , CNRS /Universite Joseph Fourier, BP 166, F-38042 
Grenoble Cedex 9, France e-mail: Bart.van-Tiggelen@grenoble.cnrs.fr 

2 Laboratoire de Physique des Plasmas, des Materiaux Conducteurs, et leurs Applications, Department of Physics, Universite 
des Sciences et de la Technologie d'Oran, USTO, Oran 31000, Algeria 

Received: date / Revised version: date 

Abstract. In this work we re-formulate and solve the self-consistent theory for localization to a Bose- 
Einstein condensate expanding in a 3D optical speckle. The long-range nature of the fluctuations in the 
potential energy, treated in the self-consistent Born approximation, make the scattering strongly velocity 
dependent, and its consequences for mobility edge and fraction of localized atoms have been investigated 
numerically. 
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1 Introduction 

Anderson localization is by now a phenomenon that has 
been widely investigated, both theoretically and exper- 
imentally, and for many different kinds of waves, from 
electrons, to electromagnetic waves, ultrasound pQ, and 
cold atoms [2]. To understand Anderson localization and 
to provide quantitative predictions for experiments, many 
tool models have been proposed, among which the Ander- 
son model is undoubtedly the best known. This model de- 
scribes a noninteracting and electron, tightly bound to the 
nucleus, but capable to tunnel to nearby atoms. Already 
in the celebrated 1958 paper [3], Anderson demonstrated 
how this model highlights the role of dimensionality. A 
genuine mobility edge only occurs in dimensions larger 
than 2 0]. For classical waves, a few observations on 3D 
localization have been reported [5l[6]. 

The tight binding model is highly relevant to under- 
stand 3D dynamical localization of cold atoms [7]. How- 
ever, is not appropriate to describe localization of many 
other waves, where the starting point is much more a dif- 
fuse, extended motion, rather than a tightly bound state. 
Localization of electromagnetic waves in 3D disordered 
media for instance, is much more a diffusion problem than 
a problem of tunneling to nearest neighbors. The scaling 
theory of localization [8] , as well as the Thouless criterion 
[9], both deal with conductance, and thus use the diffuse 
motion as a starting point. 

The self-consistent theory, first formulated by Voll- 
hardt and Wolfle in 1981 for 2D electron conductivity [10], 
was the first work that explicitly calculated how quantum 
corrections affect the classical "Drude" picture, to make 
way for localization. Despite its evident perturbational na- 



ture, the theory has been successful because it provides a 
microscopic picture of finite-size scaling, reproduces the 
Ioffe-Regel criterion for the mobility edge in 3D, and lo- 
cates the mobility edge of the tight-binding model quite 
accurately [TT] . The aim of the present work is to revisit 
and apply this theory to the localization of cold atoms in 
optical speckle. 

The first experiments on ID cold atom localization 
have been carried out recently [12] . The atoms are re- 
leased from a BEC and subsequently expand in a poten- 
tial energy landscape created by optical speckle, supposed 
free of mutual interactions. Both theory [13j and experi- 
ment have revealed the presence of a quasi-mobility edge 
in ID. Atoms with velocities v > fi/m^ (£ is the correla- 
tion length of the disorder) can hardly be scattered be- 
cause this would imply a momentum transfer larger than 
h/£ which the random speckle cannot support. As a result 
the localization length is infinite in the Born approxima- 
tion, though finite and large when all orders are taken into 
account. This somewhat surprising result highlights the 
impact of long-range correlations in ID. In higher dimen- 
sions, small angle scattering can still lead to small enough 
momentum transfer to be transferred to the speckle, even 
for large velocities, so that this quasi-mobility edge does 
not occur. Yet, correlations are expected to affect localiza- 
tion, since the potential field sensed by the atom strongly 
depends on its velocity, and strong forward scattering is 
not favorable for localization to occur. In addition, near 
the 3D mobility edge the disorder is necessarily large so 
that the spectral function of the atoms is not strongly 
peaked near energies E = p 2 /2m, neither has it a Lorentzian 
broadening. 
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2 Self-consistent Born Approximation 

In the following we consider the scattering of a nonin- 
teracting atom with energy E and momentum p from a 
disordered potential V(r). Two properties are specific for 
an optical potential. Firstly, the fluctuations 6V(r) are 
determined by the optical intensity and not by the com- 
plex field. This means that they are not Gaussian but 
rather Poisonnian. As a result, the two-point correlation 
will in principle not be sufficient to describe the full scat- 
tering statistics. Secondly, the correlation function, given 
by (SV(r)SV(r')) = [/sinc 2 (Z\r/£), is long range. Here 
U = (V) 2 depends on the average optical intensity. It is 
not difficult to see that the Born approximation breaks 
down at energies E < U/E^ [T3], with E^ = h 2 /2m^ 2 an 
important energy scale related to correlations. As usual 
we expect matter localization to occur at small energies, 
near the band edge of the spectrum. Another consequence 
of the long-range correlations is that scattering strongly 
depends on the De Broglie wave length and thus on the 
velocity of the atom. This makes it impossible to define 
a mean free path £ in the usual way, that is from the ex- 
ponential decay of the ensemble-averaged Green function 

In the following we shall cope with the second prob- 
lem. The ensemble-averaged Green function is written in 
terms of a complex self-energy as G(E, k) = [E — p 2 /2m — 
S(E,p)]~ 1 [15]. We shall apply the Self-consistent Born 
Approximation (SCBA) according to which the complex 
self-energy S(E,p) of the atom is calculated from [16] 



27(25, = £ 



U(k - k') 



E-W 



S{E, k') 



(1) 



In this equation, ^ k = J d 3 k/(27r) 3 , and from now, all 
energies, including S(E,p), are expressed in units of the 
energy scale E$. Momenta are expressed as p = fi.k with 
the De Broglie wave number k expressed in units of l/£. 
In Eq. |T]), U(k — k') represents the structure function 
associated with the speckle correlation, which determines 
the angular profile in single scattering [14J . The SCBA is 
convenient because its imaginary part expresses the gener- 
alized optical theorem in single scattering [15]. In addition, 
it avoids the bound state at negative energies predicted by 
the first Born approximation. In Ref. [T7] the SCBA was 
solved analytically for cold atoms and zero-range correla- 
tions. 

Equation fl} has been solved by iteration, with spline 
interpolation between 500 points < k n < 3. The angular 
integral can be performed analytically. Typically 10-20 it- 
erations have been necessary to ensure good convergence. 
In Figure 1 we show real and imaginary part of £(E, k) 
for U/E 2 — 1 and E — 0, and compare it to the first order 
Born approximation (FBA) applied in Ref. [14] . As a real- 
istic experimental reference we take 87 Ru-atoms released 
from a BEC with chemical potential /i = 219 h Hz into 
an optical speckle with correlation length £ = 0, 26 fim. 
This reveals that /i and E^ are equal energy scales in typ- 
ical experiments. Equivalcntly, the De Broglie wavelength 
and the correlation length are competing length scales, 



A/27T > £. To discriminate "trivial trapping " in deep ran- 
dom potential wells from "genuine" Anderson localization, 
experimentalists wish to arrange the experiment such that 
the typical kinetic energies [x are somewhat larger than 
the typical fluctuations \/U in the potential energy |12j . 
In that case, U/E 2 > 1. We will comment on this choice 
later. 
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Fig. 1. Imaginary (top) and real part (bottom) of the self 
energy (in units of E^) of an atom in a speckle potential with 
U/E^ = 1, as a function of wave number (in units of l/£) 
for an energy E = 0, calculated in the self-consistent Born 
approximation (SCBA). The dashed line denotes the first Born 
approximation (FBA) . To compensate for the shift in the band 
edge (E b /E ( = -0.85) predicted by the SCBA, the FBA has 
been evaluated at E/E^ = 0.85. 



Figure 1 (top) shows — Im ^(0, fc) to be nonzero only 
for k < 3. It is also seen that the FBA significantly overes- 
timates the amount of scattering. The real part Re S(0, k) 
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is clearly negative. This shifts the band edge of the en- 
ergy spectrum to = —0.85. The energy spectrum has a 
typical lower bound —VU, but the SCBA always locates 
the band edge at somewhat larger energies. The SCBA 
does probably not treat the (small) density of states near 
E w — \/U very well, where sharply localized Lifshits- 
type states are likely to occur. Figure 2 shows the dis- 
tribution of wave numbers at energy E = 0, expressed 
by the spectral function S(E,k) ~ —lmG(E,k). It is a 
rather broad distribution, with a tail extending to k = 2. 
This is important, since we will see in the next section 
that atom transport is quite sensitive to large momen- 
tum transfers, involving large fc-vectors. At low energies, 
relevant for localization, h/£ has become the typical mo- 
mentum of an atom, and Eg the typical energy. At larger 
energies E = 3, the spectral function behaves normally, 
i.e. strongly peaked near k = \f~E~. The peak is smaller 
because we did not conclude the geometric Ank 2 surface 
factor in phase space, so to highlight its weight at small 
k < 1 ("slow atoms") for later purposes. 




Fig. 3. Model of the atomic expansion. The BEC releases 
noninteracting atoms with kinetic energies between and its 
chemical potential fj,, that penetrate the disordered speckle po- 
tential, with average and its typical fluctuations equal to yU. 
In this work localized states (red) are found only for energies 
below the average. 
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Fig. 2. Spectral function, normalized to the total density of 
states at energy E, calculated in the SCBA, for E — (red) 
and E = 3 (blue) and for U/E^ = 1. At small energies it has 
a large weight at k — (very cold atoms) and extends up to 
fe = 2. 



3 Bethe-Salpeter equation 

We proceed with the calculation of the diffusion constant 
of the cold atoms, and the possible presence of a mobility 
edge where it vanishes. With that information we will find 
how many atoms will be localized. The idealized model we 
consider is schematically drawn in Figure 3. 

The Bethe-Salpeter equation is a rigorous equation for 
the two-particle Green function [18]. In phase space this 
object is written as <P] i ] i '(E,t,r), which is readily inter- 
preted as the "quantum probability density" for an atom 



with velocity Kk./m to travel, during the time t, from po- 
sition r = to position r and to achieve the velocity 
Kk' /m. Its Fourier-Laplace transform with space-time is 
written as ^kk' (E, f2, q). It has two fundamental prop- 
erties, namely reciprocity, <P^(E,f2,q) = < P\ l '] l (E, 12, q), 
and normalization. The last property can be expressed by, 



,_ _ n . -2lmG(E,k' 
}_^$w(E,f2,q = 0) = — ^ 



In particular, when also integrating over all energies, 



(2) 



1 



k.k 

oo 



dE—lmG{E,k) = -L- (3) 

lit /_-_ 7T — ill 



The last equality, needed for later purposes, follows from 
a general sum rule of the spectral function [18] . This iden- 
tity guarantees that the total quantum probability for 
the atom to be somewhere, with some velocity and with 
some energy, is conserved, and equal to one. The Bethe- 
Salpeter can be re-written as a quantum-kinetic equation 
for ^kk'(-E"j q)- The conservation of quantum probabil- 
ity guarantees the existence of a hydrodynamic diffusion 
pole. We shall express this as 



where 



<K£,k,q)<K£,k',q) 



^2 k -ImG(E,k) -iS7 + D(E)q 2 



(4) 



4>{E, k, q) = -Im G(E, k) - <(k • q)F(E, k)+0{q 2 ). 
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This expression states that the distribution of atoms 
in fc-space at given energy E is essentially governed by the 
spectral function, with a small correction that supports a 
current. The front factor in Eq. ((4| is imposed by the nor- 
malization condition The k- integral of —ImG(E, k) is 
recognized as (n times) the density of states (DOS) per 
unit volume p(E). With the correct normalization, we can 
set Q = 0. The still unknown function F(p) follows from 



F(E, k) = \G(E,k)\ 



2 dReG{E,k) 
dk 2 



|G(Ei ^ E »£^^ it , (E|0) 



dk' 



k 2 



kk' 



+ \G(E,k)\ 2 Y^F(E,k!)—^Uw{E^) (5) 
k' K 

The irreducible vertex J7kk' (E, q) generalizes the function 
U(k — k') defined in the first section to all interference 
contributions in multiple scattering. Once we have solved 
for F(E, k), the diffusion constant follows from the Kubo- 
Greenwood formula [TS1. 



D{E) 



h 2 l 



, 3 wp(E) 



^k 2 F{E,k) 



(6) 



Note that D(E) is determined by the fourth moment of the 
distribution F(E, k), which puts a large weight on "fast" 
atoms. The order of magnitude of the diffusion constant 
is governed by the ratio h/m of Planck's constant and the 
mass of the atom, the second factor being dimensionless 
and of order unity at low energies. For 87 Ru, h/m 1800 
pm 2 /s. The third term in Eq. (|5|) can be transformed using 
the exact Ward identity, 



E(E, k + iq) - E*(E, k - iq) = J^Uw (E, q) 



Three levels of analysis exist. First, in the Drude ap- 
proximation, one neglects all contributions form the BS- 
equation and one adopts F(E,p) — Fq(E,p), including the 
wave number derivatives. The Drude diffusion constant 
Dd{E) can be used to define a dimensionless Ioffe-Regel 
type parameter from the relation Dd{E) = ^(2h/m)k£. 
Hence, 

ht _ ^F Q (E,k) 
kl ~j: k -lmG(E,k) W 

For a short-range correlation, the self-energy is indepen- 
dent of k so that F Q (E, k) = 21m 2 G(E, k), and this defini- 
tion of kl coincides with the usual one in terms of Im S(E) 
|15j . For low energies we found in Figure 3 that typically 
few l/£ . If we anticipate that ki ~ 1 near the mobility 
edge, we conclude that the mean free path is roughly equal 
to the correlation length. 
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x ^G(E,k' + ^ q )-G*(E,k'-^ (7) 

If this identity is developed linearly in q, and inserted into 
Eq. ([5]) , we obtain, 



F(E, k) = F (E, k) + SqU{E, k) 

+ \G(E,k)\ 2 J2nE,k')^U kk >(E,0) (8) 



with 



F (E,k)= \G(E,k)\ 2 1 + 



o/, , dReE{E,k)\ 8ReG{E,k) 



dk 2 



\G{E,k)\ 2 ^lmG{E,k') hi 



' k 2 



d_ 
dq 



Ok 2 



Fig. 4. Solution of the BS equations ((6]) and ((8| in the Boltz- 
mann approximation, for a disorder amplitude U/E^ = 1. 
Shown as a function of energy are the Drude conductivity 
Tvp(E)D d (E), the Ioffe-Regel parameter |9j|, the DOS that van- 
ishes at the band edge E/E^ = —0.83, and the ratio of Boltz- 
mann and Drude conductivity. Note that 2/3 times the Ioffe- 
Regel parameter equals the Drude diffusion constant, expressed 
in units of h/m. 



The Drude approximation is popular in electron - im- 
purity scattering but clearly inadequate when the scatter- 
ing is strongly anisotropic, as for the optical speckle. In the 
Boltzmann approximation we adopt C/kk' = U(k— k ), i.e. 
the structure function associated with the optical disorder. 
Being a function of k— k' only, it follows that 5qU = 0. Dif- 
ferent results are summarized in Figure 4. At low energies 
Boltzmann and Drude diffusion constant typically differ 
by a factor 1.5 as was obtained by Kuhn etal [2] on the ba- 
sis of the FBA. For atom energies E > E% the Boltzmann 
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0,5 1,0 1,5 2,0 2,5 3,0 

Wave number k 

Fig. 5. The function fc 4 Fs(_E,fc), solution of the BS-equation 
for the energy E/E^ = —0.06 in the Boltzmann approximation. 
The blue dashed line compares it to the Drude approximation 
k 4 Fq(E , p) defined in Eq. JS), rescaled properly such that both 
k- integrals equal as/cro 



diffusion constant rapidly rises since only strong forward 
scattering can occur. In the region E = 0, the Ioffe-Regel 
parameter takes values of the order of 1.5. We infer from 
Figure 5 that the solution Fb(E,p) is roughly a rescaling 
of the Drude Ansatz, for which the current is dominated 
by atoms with velocities v — 0.83?i/m£. Nevertheless, a 
non- negligible fraction of atoms faster than v = l.bTi/m^ 
contributes to the current (note that fi/m^ « 7 mm/s for 
the set-up with Rubidium described above). 

We finally consider constructive interferences, and add 
the most-crossed diagrams to the BS-equation in the spirit 
of the self-consistent theory of localization. Any observa- 
tion D < Db must be attributed to constructive interfer- 
ences. It is well-known that, by reciprocity, these diagrams 
can be constructed from the solution 4>(E, k, q) of the BS- 
equation by removing the incoming and outgoing Green's 
functions (indicated by a hat), and by time-reversing the 



bottom line [19] 



k'). Using Eq. (j4|) the two procedures lead to 



#kk'(£) 



a{EW 

[ImE(E,k)lmE{E,k') + tq ■ {hk< {E)k + I k , k (E)k')} 
U^(E, q} 



a(E)(k + k>) 



1 



Im 2 ^, -Ak) + I A k/2,Ak/2{E)i(k + k') ■ q 



(10) 



We abbreviated I kk > = F{E,k)ImS(E,k')/\G(E,k)\ 2 , Ak 
= |k— k'| and introduced a{E) = irp(E)D(E), the equiv- 
alent of DC-conductivity in electron conduction. We now 
face the more complicated task of solving Eqs. ((8]) and 



© simultaneously with J7 kk / = U(k - k') + and of 

finding out if its extrapolation to small energies leads to 
a mobility edge where a(E) = 0. This constitutes a "self- 
consistent" problem for the entire function F(E, k), rather 
than just for its fourth moment, the DC-conductivity. 

We first observe that for most-crossed diagrams SqU ^ 
0. This term does not appear in standard self-consistent 
theory [10] . which relies on moment expansion. Note that 
it also features " self-consistently" the function F(k)/a, 
just like the second term in Eq. ((8]). As can be induced 
from Eq. (fTU|) the singularity at k = — k' that generates 
the (weak) localization is partly compensated by the factor 
k + k'. We shall therefore ignore it here as well. 

The self-consistent equations (jU),© and (p~0|) can be 
solved almost analytically when the self-energy is assumed 
fc-indepcndent, typically true for zero-range correlations. 
Without more details we mention that the mobility edge 
then occurs at ki — \.122 [20 . Quite convenient is that, 
even when scattering extends to infinite wave numbers, 
our theory does not require an ad-hoc cut-off to elimi- 
nate short wave paths that diverge in approximate theories 
[CTET] . For the speckle correlation, Figure 6 gives the re- 
sult of the exact numerical solution, obtained by iteration 
and spline interpolation. This method worked satisfacto- 
rily until close to (AE w E^/5) the mobility edge, where 
the most-crossed diagrams give a diverging contribution. 
Before that happens, the strong forward scattering of a 
single scattering competes heavily with the reduction in- 
duced by weak localization. 



U&?(E, q) = % k -k'+ q) i(k'-k+ q )(£, k+ 




-0,8 -0,6 -0,4 -0,2 0,0 0,2 0,4 0,6 0,8 1,0 

Energy / 



Fig. 6. Solution of the BS equation (|8]) with inclusion of the 
most-crossed diagrams for a disorder amplitude U/E^ — 1. 
Shown in red is the ratio of conductivity and Drude conductiv- 
ity for energies E > 0.15-B^ for which our iteration converged. 
The blue dashed line relies on an approximation discussed in 
the text and locates the mobility at E = — 0.06-E^. 



To find the location of the mobility edge we shall use 
the following approximation. In Figure 7 we see that at 
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Fig. 7. The function k 4 F(E,k) (in red)), solution of the BS- 
equation for the energy E/E^ = 0.15, compared to the Drude 
approximation k 4 Fo(E,p) defined in Eq. (|5J (blue dashed line). 
The first has been rescaled by the factor ag/a such that the 
fc-integrals of both are equal to one. 

E/E^ = 0.15 the solution F(E,k)/a is closely approx- 
imated by Fo(E,k)/cro (both have their fourth moment 
normalized to one) . This equivalence is physically reason- 
able since it implies that all atoms with energy E undergo 
the same reduction in diffusion, but keep the same veloc- 
ity distribution as found in the Drude picture. If we insert 
Fq (E, k) in the left hand side of Eq. and integrate over 
k we can derive the simple relation <j/<7b ~ 1 — K{E) that 
is reminiscent of standard self-consistent theory [TO] ■ The 
parameter K is found to be, 

Im 2 i;(£,i|k-k'|) \G(E,k)\ 2 (11) 

For U/Ec = 1, this approximation locates the mobility 
edge (K = 1) at E/E^ = -0.06 (0.94 when we shift ener- 
gies over \JT] as in Fig 8). Upon inspecting the numerically 
exact solution in Figure 8, we suspect that the real mo- 
bility edge is located somewhat lower, near E/E^ = —0.1. 
For smaller disorder we found that the approximation be- 
comes better. It is in principle possible to obtain numeri- 
cally the function F(E,p)/a as a — > and calculate more 
precisely the location of mobility edge, but this is beyond 
the scope of this paper. In Figure 8 we show the different 
energies for different strengths of the disorder. If we apply 
the criterion found in Ref [13) . Mb = 0.95 (dashed line 
in Figure 8), localization would occur at smaller energies, 
around E/E 6 = -0.35 (0.65 in Fig. 8) for U/E% = 1. This 
approach expresses the general trend well but is clearly 
somewhat pessimistic, probably because their choice of 
the ad-hoc wave number cut-off to calculate the most- 
crossed diagrams underestimates K. Note that all local- 




Disorder U / E 



Fig. 8. Band edge and mobility edge as a function of the disor- 
der strength U, all expressed in the energy E^. Contrary to the 
former figures the energy scale has been shifted by the average 
potential yfU of the potential energy so that a direct compar- 
ison can be made to the chemical potential > 0, and the 
criterion E c — E b w 0.5U/E^ found by Ref.[l4] (blue dashed). 
The numbers around the mobility edge reveal the small vari- 
ation of the Ioffe-Regel parameter ([9|. The exclusion of small 
energies is likely to be an artifact of the SCBA and we expect 
strongly localized Lifshits-type states with small DOS. 

ized states occur below the average of the potential land- 
scape. Unlike in ID, we find no regime with atoms fast 
enough to traverse the potential barriers, but to become 
localized purely by constructive interferences, without the 
assistance of tunneling. 

A final important question is how many atoms will 
be localized, given an initial velocity distribution that is 
determined by the expansion of the BEC after eliminat- 
ing the trap. We emphasize that according to the scenario 
sketched in Fig. 3, the chemical potential [i of the BEC 
does not represent the energy of the atom inside the dis- 
ordered potential, but rather the distribution of incident 
kinetic energies. If this distribution is denoted by <^> M (k), it 
follows from Eq. ((4]) that the fraction of localized atoms, 
regardless of their final velocity or position, is given by 




= f ' gIEY^— ImG(£,fc> M (k') (12) 
Je„ k , t 

Castin and Dum [52] showed that after the free ex- 
pansion 4>fj,{k) ~ 1 — k 2 /k^ and zero for k > fc M , where 
the maximum wave number k^ — WJl/E^. The fraction 
of localized atoms is thus determined by the number of 
microscopic states below the mobility edge whose kinetic 
energies are smaller than /x. In this discussion it is conve- 
nient to make the zero of the energy scale the same for /i 
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and E, as was already done already in Figure 8, shifting 
the localized region to positive energies. In present (ID) 
experiments is U/E^ = 1 and /i « Eg which is only slightly 
above the 3D mobility edge. One might perhaps expect 
most atoms to be localized. Our calculations clearly show 
that the distribution of atom energies is quite different 
from ^(k). This (normalized) distribution F(E) is given 
by the wave number integral in Eq. (|10p and is shown in 
Figure 9. For /i <C Eg it reduces to the spectral function 
at k = 0, independent of M (fc). Even for small /i, many 
atoms achieve energies E > E c (40% for U/E^ = 1) and 
are delocalized. This number agrees well with predictions 
based on zero-range correlations in which case 45 % was 
found to be localized for /i <C E c [17] . The fraction of delo- 
calized atoms further decreases as the chemical potential 
rises (Figure 10), with only 35% localized for fj, = Eg. 

Even if we choose /i > ^/U, the atoms that localize have 
energies E < s/XJ . 



0,6- i\ 




chemical potential u/E^ 

Fig. 10. Fraction of localized atoms as a function of chemical 
potential, for different amplitudes of the disorder. 



1,0 
0,9 
0,8 

C 

O 0,7- 



~& 0,4 

0,3 

<D 

C 0,2 

Q) 

0,1 
0,0 



"H 1 



I \ \i = 0.25 E E 

f : \ V 
■ i .»•••■•••• 



;/ i <^ Xr 2E * 



0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 



energy /E^ 



Fig. 9. Energy distributions of the atoms inside the disordered 
speckle potential, for different chemical potentials of the BEC 
from which they were released and U /E^ = 1. They all exhibit 
a tail of relatively fast atoms (E > E$) that extends beyond 
the mobility edge E c . 



In conclusion, we have calculated the phase diagram 
for localization of cold atoms in a 3D speckle potential, 
using the self-consistent Born approximation and the self- 
consistent theory of localization. The mobility edge is char- 
acterized by a Ioffe-Regel type parameter that varies be- 
tween 1.2 and 1.4. Depending on the chemical potential 
of the BEC, typically 35% to 60% of the atoms are lo- 
calized. The self-consistent Born approximation deals al- 
ready much better with the long-range correlations and 
the broadness of the spectral function than the first Born 
approximation, but does not discriminate between differ- 
ent statistics of the disorder. Yet, the mobility edge of 
the tight-binding model is known to depend on that [11] . 
It would be very interesting to apply a recently proposed 
method [53] to calculate the self-energy of the atoms more 



precisely. The theory presented in this work can then be 
used straightforwardly to find the mobility edge. 

We would like to thank Sergey Skipetrov for help and discus- 
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